Nonlinear optical response and spin-charge separation in one-dimensional Mott 

insulators 
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We theoretically study the nonlinear optical response and photoexcited states of the Mott insulators. 
The nonlinear optical susceptibility x^^^ is calculated by using the exact diagonalization technique on 
small clusters. From the systematic study of the dependence of x'^' on dimensionality, we find that 
the spin-charge separation plays a crucial role in enhancing x*"^' in the one- dimensional (ID) Mott 
insulators. Based on this result, we propose a holon-doublon model, which describes the nonlinear 
response in the ID Mott insulators. These findings show that the spin-charge separation will become 
a key concept of optoelectronic devices. 

PACS numbers; 78.20.Bh, 42.65.-k, 71.10.Fd 
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The charge gap in Mott insulators is a consequence of 
strong electron correlation. This is completely different 
from the band insulators, where the charge gap is basi- 
cally originated from band effects. The nature of charge 
excitation across the gap is thus essentially different be- 
tween the two types of insulators. The optical response 
is used to investigate the charge excitation across the 
gap. In the response, the linear susceptibility with re- 
spect to the applied electric filed, x^^"*; provides infor- 
mation on the dipole-allowed states with odd parity. In 
addition to x^^^ the nonlinear susceptibilities detect the 
odd states together with the dipole-forbidden states with 
even parity .cl 

Very recently, anomalously enhanced third-order non- 
linear optical susceptibility x'^' has been reported for 
one-dimenajonal (ID) Mott insulators of Cu oxides and 
Ni halides,a as compared with those for the band insula- 
tors. In addition, the ID Cu oxide, Sr2Cu03, exhibits 
ultrafast nonlinear optical response (~1 ps) at room 
temperatureE These facts strongly suggest a great po- 
tential of the ID Mott insulators as optoelectronic mate- 
rials with high performance J3'u Furthermore, it is shown 
that x*^'^-' in 2D Mott iiasulators is smaller than that in 
the ID Mott insulators J3 

Motivated by these striking experiments, we theoreti- 
cally examine photoexcited states and nonlinear optical 
response in the Mott insulators, and clarify underlying 
physics of optical nonlinearity of the ID Mott insula- 
tors. We use the half- filled Hubbard model to describe 
the Mott insulators. The x^^^ is calculated by using the 
numerically exact diagonalization method on finite-size 
clusters, and compared with that for higher dimensional 
systems with ladder and 2D geometry. We find that 
X^'^^ increases with decreasing the dimensionality. It is 
shown that in the ID system dipole-allowed (odd) and 
-forbidden (even) states are almost degenerate in energy, 
having very large dipole coupling between the two states, 
while in the 2D system the dipole coupling is rather small 
in spite of the closeness of the odd and even states. We 
demonstrate that this peculiar feature in the ID system 
is due to the spin-charge separation inherent in the ID 



Mott insulators. We propose an effective model that can 
describe the optical nonlinearity, a holon-doublon model, 
where holon and douhlon represent the charge degree of 
freedom for photoinduced unoccupied and doubly occu- 
pied sites, respectively. The model reproduces very 
the characteristic behaviors of the experimental x'''^'' 
We find that the spin-charge separation becomes a key 
concept of future optoelectronic devices. 

The electric field E of the incident light induces the di- 
electric polarization Pind in a material, which is described 
by a power series of nonlinear optical susceptibility x*-"-*: 
Pind = eo (x^^^E + x^^^E^ + x^^^E^ 4- • • •). The linear 
susceptibility x*-^-* is given by 



E 



'Oa'aO 



'Qa' aQ 



ria + iTa 



(1) 



where L is the number of sites, eq is the dielectric con- 
stant, j and k are the polarization directions, er^a is the 
dipole moment between the ground state |0) and excited 
state I a) with odd parity, Qa is the energy difference be- 
tween |0) and |a), and Ta is the damping factor. Due to 
symmetry restrictions, x'^' vanishes in centrosymmetric 
materials to which Cu oxides and Ni halides belong. The 
lowest observable nonlinearity is, therefore, x^"^^) which 
is expressed ai 
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where lo^ = uji -\- UJ2 + i^s, b and c denote even and 
odd states, respectively, and V represents the sum of 
permutation on (j, wi), (fc, 102), {I, ^^3), and (m, 
LUa-). In the present study, by setting all of the polar- 
ization directions the same along chains, we examine 
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parts of which give two-photon absorption (TPA) and 
electroabsorption spectra, respectively. Hereafter, we set 
e, h, and eo to be unity. 

The insulating Cu oxides and Ni halides are known 
to be charge-transfer (CT)-type Mott insulators, where 
both 3d and 2p orbitals exist in the transition metal and 
ligand ions, respectively, and participate in the electronic 
states. The values of the gap are predominantly deter- 
mined by the energy position of the p orbitals. However, 
it is well established that the electronic states of the CT- 
type insulators can be described by the Hubbard model 
with single band by mapping a bound state called the 
Zhang-Rice singlet state onto the lower Hubbard band. 
Actually, the single-band Hubbard model explains very 
well the spectral line shape of angle-resolved photoemis- 
sion and electron^egergy loss spectroscopies in the insu- 
lating Cu-oxides.ETtl 

The single-band Hubbard Hamiltonian is defined as 

H = -t (4„Cj-^ + H.c.) +;7^n,,|n,,i 
+ V"^n^nj, (3) 

where cj ^ is the creation operator of an electron with spin 
a at site i, ni—ni^-i'+rii^i, runs over pairs on nearest 
neighbor sites, t is the hopping integral, U is the on-site 
Coulomb interaction, and V is the Coulomb interaction 
between nearest-neighbor sites. |-| r-. 

The realistic values of the parameters for cupratestlcl 
are chosen: U/t=10, and V/t=1.5 for a ID system and 
V/t=l for ladder and 2D systems. Additionally, in 
the 2D system, a second-nearest-neighbor hopping t' is 
considered to simulate a parent compound of high- Tc 
cupratesB'B The damping factors are assumed to be tlie 
same for all excited states of the systems with r=0.4t.EI 

We numerically examine linear absorption spectra, 
that are defined by an imaginary part of x^'^\ ^^'^ TPA 
spectra for a 12-site chain with open boundary condition. 
The Lanczos technique is used for the calculation of the 
ground state. The spectraj-are calculated by using the 
correction vector technique.EB We also calculate them for 
two-leg ladder, and 2D systems with 12 sites in order to 
clarify the dependence of x'^' on dimensionality. Thfi. 
clusters used in the calculation are shown in Fig. 1(c) .Ell 

The calculated results are shown in Figs. 1(a) and 1(b). 
The solid, dashed, and dotted lines denote the results for 
the ID, ladder and 2D systems, respectively. The upper 
panel (a) represents the linear absorption spectra, which 
detect odd states. The optical gap is estimated to ^2.4 
{—6t) eV assuming 1f4^A eV, comparable with the exper- 
imental optical gapJl3 The linear absorption has a large 
magnitude for the ID system with strong enhancement 
in the low-energy region {ui ~ 7t), and decreases with 
increasing dimensionality from the ID to 2D systems. 
These results are consistent with the experimental data 
of optical conductivity for ID, ladder, and 2D cuprates.Ej 
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FIG. 1. The linear absorption (a) and TPA (b) spectra for 
the ID, ladder, and 2D Hubbard models, denoted by solid, 
dashed, and dotted lines, respectively. The clusters used in 
the calculation are shown in (c). U/t— 10. V/t=l.b for the 
ID cluster and V/t=l for the ladder and 2D clusters. For the 
2D cluster, t' /t^-OA. 

The lower panel (b) shows the TPA spectra, which 
detect even states with the same parity as the ground 
state. We find that the spectral shapes of TPA bear 
a close resemblance to those of x^^^ when the energy 
in TPA is doubled. This means that the energy posi- 
tions of the even states lie in the same energy region as 
the odd sta|tes shown in Fig. 1(a) [see the denominator 
of Eq. (2)].E£l These characteristics are independent of 
the dimensionality, and agree with the results obtained 
by pump-probe measurements. B'Q We also find that the 
magnitude of x^'^^ is the largest in the ID system, and de- 
creases with increasing the dimensionality. This indicates 
that one dimensionality is favorable to the enhancement 
of x^^^ I which is consistent with th&|recent experimental 
data on Sr2Cu03 and Sr2Cu02Cl2.Q 

In order to clarify the TPA spectra in detail, we next 
examine photoexcited states and dipole moments be- 
tween these states. The dipole moments are given by 
^j(odd|a;irti|even), where |odd) and jeven) denote eigen- 
states with odd and even parities for the Hamiltonian of 
Eq. (3), respectively, Xi is the x-coordinate of the posi- 
tion at the i site, and rii is the number operator of the 
electron. 

Figure 2 shows the results of the energy distribution 
of dipole moments for the ID (left panel) and 2D (right 
panel) systems, which are calculated exactly by using 
clusters shown in the insets. The upper panel [(a) and 
(d)] represents the dipole moments between the ground 
state and odd states. The middle panel [(b) and (e)] 
shows dipole moments between the first odd state labeled 
as Oi in (a) and (d) and even states. In the ID case, the 
Oi state strongly couples to an leym state with nearly 
the same energy as the Oi stateO'til The magnitude of 
the dipole moment is about 3. On the other hand, in 
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FIG. 2. The energy distribution of the dipole moments for 
8-site ID [(a), (b), and (c)] and 2D clusters [(d), (e), and (f)]. 
The (a) and (d) show the dipole moments between the ground 
state and the dipole-allowed states. The (b) and (e) [(c) and 
(f)] show the dipole moments between the dipole-forbidden 
states and typical two dipole-allowed states, denoted by 0\ 
(O2) shown in (a) and (d). The S denotes a quantity to mea- 
sure the similarity of the wave functions. )7/t=10. l^/t=1.5 
for the ID cluster, and F/t=l and t' /t=-OA for the 2D clus- 
ter. 



the 2D case the dipole moment is smaller than that in 
the ID system and the magnitude is less than 0.5. In 
addition, other dipole moments are also very small in the 
2D system. The lower panel [(c) and (f)] shows dipole 
moments between another odd state O2 and even states. 
We find again that an even state with very close energy 
to the O2 state yields large dipole moment in the ID 
case, whereas the dipole moment in the 2D case is small 
although there are even states whose energy is very close 
to the O2 state. These results indicate a quantitative 
difference of the dipole moments between the ID and 2D 
systems. The difference plays an important role in the 
difference of the magnitude of TPA between the ID and 
2D systems. 

In general, the dipole moment becomes large when two 
states which lead to the dipole moment are constructed 
by similar wave functions. Therefore it is meaningful to 
examine the nature of the odd and even states in the 
ID and 2D systems. We define a quantity S, which 
measures the similarity of the wave functions, as follows: 
S = |ao,iae,i|, where Uo^i (ae.i) is the coefficient of a 
basis which constructs the odd (even) eigenstates of the 
Hamiltonian of Eq. (3), and i is the index of the bases. 
The quantity S becomes unity when two eigenstates are 
the same except for the phase of each basis. 

In the ID system, S between the Oi (O2) state and an 



adjacent even state is 0.96 (0.81), which is close to unity. 
On the other hand, in the 2D system, S is 0.66 and 0.64 
for the Oi and the O2 states, respectively. These results 
show that in the ID system the two wave functions with 
close energies are quite similar to each other in contrast 
to the 2D ones. 

Next, we consider the origin of the different behavior of 
the photoexcited states between the ID and 2D systems. 
In the ID system, a photoinduced carrier separates into 
two collective modes carrying the spin and charge de- 
grees of freedom (spin-charge separation). In addition, 
the wave function itself is factorizcd as the products oL 
spin and charge wave functions when U is-^ery large .Ej 
Therefore, it is reasonable in the ID systemll3 to take into 
account only the charge degree of freedom of electrons 
for the optical response. We propose below an effective 
model to describe the nonlinear optical response in the 
ID system. In contrast to the ID system, the motion 
of the carriers in the 2D system is known to be strongly 
influenced bypSpin background as long as the excitation 
energy is 1ow,Lj implying that the wave function cannot 
be simply factorizcd unlike the ID system. It is, there- 
fore, natural to suppose that the spin degree of freedom 
in the 2D system plays a crucial role in the broad spec- 
tral weight of x^^^ and x^'^'' (Fig- 1) as well as the broad 
distribution of the photoexcited states (Fig. 2). 

Based on the spin-charge separation picture, we intro- 
duce a holon- doublon model for the ID Mott insulators. 
The Hamiltonian is given by 

H = tJ2 ~ dld,+i + H.c.) -VY^ h\hd]dj 

■Y,{h\h+4d.), (4) 
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with the constraint of no double occupation at site i, 
{h\hi + d\di) < 1. The holon h (doublon d) represents 
the charge degree of freedom of unoccupied (doubly oc- 
cupied) sites produced by the photoexcitation. We note 
that only two particles, i.e., one holon and one doublon, 
are contained in the present system. Although the model 
includes the attractive Coulomb interaction, we note that 
the model is different from a model for standard exciton 
of ID semiconductors, since the holon and the doublon 
cannot occupy the same site due to the hard core con- 
straint unlike the semiconductors. As will be shown be- 
low, the constraint plays an important role in the degen- 
eracy of the even and odd photoexcited states. 

It is easily shown in the holon-doublon model that 
the even and odd states are degenerate and thus the 
overlap integral S between them is unity: Let us in- 
troduce a set of the bases with odd and even pari- 
ties written as | basis, o) = ( | hdQQ . . .) — 1 00 . . . dh) ) j \/2 and 
[basis, e) = (|/idOO. . .) -f- \00 . . .dh))/V2, where h, d and 
denote holon, doublon, and vacant sites, respectively. 
Since there is no exchange between holon and doublon 
due to the hard core constraint, the first terms of the 
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FIG. 3. The linear absorption (a), TPA (b), and elec- 
troabsorption (c) spectra obtained by holon-doublon model 
for an 80-site chain with open boundary condition. U /t=10, 
V/t=1.5, and r/t=OA. 



bases never couple to the second terms. This means that 
the matrix elements of the Hamiltonian of Eq. (4) are 
independent of the sign of the second terms. There- 
fore, the Hamiltonian matrix with odd parity is equal 
to that with even parity. This gives rise to the degen- 
eracy between odd and even eigenstates and thus 5=1. 
This model naturally explains the characteristics in the 
photoexcited states in the ID Mott insulators discussed 
above. We note that the even and odd states are always 
degenerate regardless of the value of V (Ref. 18) or the 
range of Coulomb interactions. 

Figure 3 exhibits the linear and nonlinear absorptioH. 
spectra in the holon-doublon model for an 80-site chain.tSi 
The linear absorption (a) and TPA (b) spectra reproduce 
well those for the ID Hubbard model shown in Fig. 1. 
The obtained results are also very similar to those of the 
experimental data in Sr2Cu03 (Refs. 2 and 3) in the fol- 
lowing points: (i) The spectral weights of linear absorp- 
tion and TPA spectra are concentrated on the same en- 
ergy region when the energy in TPA is doubled, (ii) The 
electroabsorption spectra (c) show an oscillating struc- 
ture at the spectral edge of the linear absorption 6t) . 
The agreement suggests that this model is a proper one 
in describing the nonlinear optical response in the ID 
Mott insulators. 

Finally we discuss the difference of nonlinear optical 
response between ID band insulators and Mott insula- 
tors. The photoexcitation process in the band insulators 
creates an electron-hole pair bound by their attractive 
Coulomb interaction (Mott-Wannier exciton). The pair 
gives a series of bound states with odd and even parities 
below the gap. IiLfact, in the ID hand insulators such as 
silicon polymersEHl and Pt halides,Ej it has been observed 



that the odd and even states are concentrated on differ- 
ent energy regions, and their difference is rather large 
(~1 eV). However, such an energy difference between 
odd and even states is generally unfavorable for obtain- 
ing the large dipole moments because dipole moments 
are inversely proportional to their energy difference. In 
this respect, the ID Mott insulators are advantageous to 
large dipole moments since they have the nearly degener- 
ate odd and even states as shown in Figs. 1 and 3, which 
leads to the large x^^^ ■ 

In summary, we have clarified the nonlinear optical re- 
sponse and photoexcited states of Mott insulators. We 
found that x*^"^^ is enhanced with decreasing the dimen- 
sionality, which is consistent with the recent experiments. 
In the ID system, the odd and even photoexcited states 
are nearly degenerate in energy, and constructed by sim- 
ilar wave functions. We found that these properties, 
which become important for the large nonlinearity, are 
originated by the spin-charge separation. We also pro- 
posed the effective model in ID Mott insulators. These 
findings imply that the concept of the spin-charge sep- 
aration is not of purely academic interest but will show 
up as underlying physics of optoelectronic devices in the 
near future. 
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